155М 1561 - 5359. Штучний 1нтелект, 2015, № 1-2 


УДК 004.93 


Л.Г. Ахметшина, А.А. Егоров 
Днепропетровский национальный университет имени Олеся Гончара, Украина 
Украина, 49050, г. Днепропетровск, ул. Научная, 13 


МЕТОД АДАПТИВНОЙ ВИЗУАЛИЗАЦИИ ПРИ СЕГМЕНТАЦИИ 
ИЗОБРАЖЕНИЙ НА ОСНОВЕ СИНГУЛЯРНОГО РАЗЛОЖЕНИЯ 


Г.@. Асйтебзйта, А.А. Уегогоу 
Ощергореноу$К Майопа! ОшуегзИу патед Бу О]е5 НосВаг, ОКташте 
ОКгатше, 49050, с. ОшергорейоузК, МаисВпаа $, 13 


ТНЕ МЕТНОР ОЕ ТНЕ АБАРТТУЕ УТЗОАТГЛААТТОМ ОЕТМАСЕЗ 
ЗЕСМЕМТАТТОМ ВЕЗОГТ$ ВАЗЕО ОМ ТНЕ УТМСОГАВ УАГОЕЗ 
РЕСОМРОЗГТТОМ 


Л.Г. Ахметшина, А.О. Егоров 
Дншропетровський нащшональний ун1верситет 1мен1 Олеся Гончара, Украйна 
Укрална, 49050, м. Дншропетровськ, вул. Наукова, 13 
МЕТОД АДАПТИВНОГ В1ЗУАЛ!ЗАЦИ ПРИ СЕГМЕНТАЦИ 
ЗОБРАЖЕНЬ НА ОСНОВТ СИНГУЛЯРНОГО РОЗКЛАДУ 


Работа посвящена описанию метода адаптивной визуализации результатов сегментации изображений, который 
основан на сингулярном разложении. Предложенный алгоритм за счет автоматизации отбора наиболее 
значимых столбцов матрицы левых сингулярных векторов и композиции с исходным снимком позволяет 
повысить наглядность и чувствительность сегментации. Представлены результаты экспериментальных 
исследований возможностей предложенного метода на примере визуализации результатов сегментации 
полутоновых и цветных медицинских изображений различной физической природы. 

Ключевые слова: обработка изображений, сегментация, нечеткая (нейро-фаззи) 
кластеризация, визуализация результатов сегментации, сингулярное разложение. 


ТЬ$ агисе 15 деуоще4 ю 4езсирноп оЁ пе тефо4 о Фе адарйуе у15иа|таНоп о{ ипасез зезтегайоп гези $ 
Базе оп Фе зшешаг уаез десотроз@оп. ТЬе ргорозе а1еогИтл аПо\уз епрапсше о @е с|еагпез$ апа 
зепзШуйу оЁ зестещаноп о\лие ю аиютанНоп оЁ Фе сБо1се оЁРШе поз $1ет1Нсап соаили$ оЁ Фе 1ей зшещаг 
уесютз таблх ап Фе тегоше Уф Фе ша! пласе. ТВе ехрегипета1 гезеагсбез оЁР фе ргорозе а]еогит 
аНиез оп ехатр/е оЁ Фе у15аНиште оЁ Фе зестешайоп гезиН$ оЁ Фе отаузсае ап4 со]ог те@са] ппазез 
уапоцз рБуз1са] пайлге аге зВо\уп. 

Кеууог@5$: пласе ргосеззте, зестетайоп, Ра7ту (пеиго-Ё12ту) стазепие, Фе у1зиаПтаНоп оЁ 
Фе зестешайоп гези $, зпеаг уаез Чесотроз!оп. 

Робота присвячена опису методу адаптивно! взуал1заци результалив сегментацйи зображень, який 
заснований на сингулярному розкладенн!. Запропонований алгоритм за рахунок автоматизаци в1дбору 
найб1льш значущих стовбщв матриц! л1вих сингулярних вектор1в та композиций з вадним зн!мком 
дозволяе шдвищити наочнсть та чутлив!сть сегментаци. Наведено результати експериментальних 
досллджень можливостей запропонованого методу на приклад! взуаллзацй результатв сегментаци 
нашвтонових та кольорових медичних зображень р1зно! ф1зично! природи. 

Ключов! слова: обробка зображень, сегментащя, нечтка (нейро-фазз1) кластеризацщя, 
в1зуаллзацля результат1в сегментаци, сингулярне розкладання. 


Введение 


Одним из наиболее сложных этапов при решении прикладных интел- 
лектуальных задач, связанных с обработкой изображений, в том числе и много- 
мерных, является их сегментация. Сложность этой процедуры обусловлена как 
существенным различием характеристик анализируемых изображений, так и 
разнообразием поставленных целей. Кроме того, отсутствие объективных критериев 
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оценок качества исходных данных и конечного результата [1] приводит к сложности 
выбора алгоритмов для выполнения этой задачи. 

В настоящее время достаточно эффективным подходом для сегментации 
изображений является выполнение кластеризации с последующей визуализацией ее 
результатов. При использовании такого способа получаемый результат существен- 
ным образом зависит как от применяемого метода кластеризации, так и способа 
формирования конечного изображения. 

Среди популярных на сегодняшний день подходов для решения задачи 
кластеризации следует выделить нейронные сети, нечеткие и нейро-фаззи 
алгоритмы. При этом особенный интерес представляют именно нейро-фаззи 
алгоритмы, т.к. они позволяют объединить достоинства нейронной сети и нечеткой 
логики в рамках одного метода. 

Многомерность используемой в нейро-фаззи алгоритмах нечеткой функции 
принадлежности, отсутствие универсального метода определения значимости 
отдельных кластеров и неоднозначность выполнения процедуры дефаззификации, а, 
соответственно, и визуализации [2] являются ключевыми проблемами, решение 
которых определяет информативность результата сегментации. 


Постановка задачи 


В работе [3] был предложен алгоритм визуализации результатов нечеткой 
кластеризации на основе применения сингулярного разложения, что позволяет 
устранять избыточность и повышать информативность исходных данных за счет 
приведения вещественной матрицы к каноническому виду, отражающему ее 
геометрическую структуру [4]. 

Суть алгоритма заключалась в том, что после осуществления сингулярного 
разложения для полученной нечеткой функции принадлежности производилась 
автоматическая оценка информативности столбцов матрицы левых сингулярных 
векторов на основе матрицы правых сингулярных векторов. Затем осуществлялся 
выбор 3-ех наиболее информативных столбцов из матрицы левых сингулярных 
векторов, и на их основе формировалось выходное изображение. 

В данной работе предложен метод адаптивной визуализации результатов 
сегментации на основе сингулярного разложения, который за счет автоматизации 
отбора наиболее значимых столбцов матрицы левых сингулярных векторов и 
композиции с исходным изображением позволяет повысить наглядность и 
чувствительность сегментации. 


Решение задачи 


Предложенный алгоритм адаптивной визуализации состоит из следующих 
шагов: 

1. Осуществляется нечеткая кластеризация исходного многомерного изображения 
[Г размерности а’хахх (4 и Ах - размерность снимка по вертикали и 
горизонтали, соответственно, (4 — количество цветовых/спектральных 
компонент) в результате чего формируется нечеткая функция принадлежности 
О размерности пхс (п=ау»-ах, с - количество нечетких кластеров). Основной 
целью этого шага является расширение пространства входных признаков 
(возможны и другие способы), что, с одной стороны, позволяет учитывать 


© Л.Г. Ахметшина, А.А. Егоров 37 


38 


Т55М 1561 - 5359. Штучний 1нтелект, 2015, № 1-2 


дополнительную информацию при сегментации, а с другой — применять 
сингулярное разложение для полутоновых одномерных изображений. 

При формировании полутонового результирующего изображения далее 
выполняются следующие действия: 


а) 


6) 


в) 


г) 


д) 
е) 


ж) 


3) 


и) 


к) 


Формируется 3-мерная матрица / ви размерностью пх(@+с), которая 
содержит исходное изображение и матрицу И. 

Полученная матрица / } подвергается сингулярному разложению [5], в 
результате чего формируются три матрицы: /' (содержит левые сингулярные 
вектора), 5; (содержит собственные числа) и Г,(содержит правые 


сингулярные вектора) такие, что: 
и=и,-5, У. (1) 
При этом размерности матриц / а р совпадают, а размерность матриц 5; и 


Г; имеет значение (+ с)х (@&+с) (используется экономное представление 


сингулярного разложения). Данный шаг обеспечивает переход к новому 
пространству информативных признаков. 


Все столбцы матрицы Г ) масштабируются на отрезок [01] и 
интерпретируются как матрица размерности 4х ах х (а + с). 

На основе матрицы Г; вычисляется вектор коэффициентов С по следующей 
формуле [3]: 


@+с @+с 
с=|| У, + > ©,), 1/2 ес], (2) 
= = 


Вектор коэффициентов С упорядочивается по убыванию, а его элементы 
нормируется так, чтобы их сумма была равна 1. 

Формируется вектор @&С , содержащий разности для каждой пары соседних 
элементов вектора С. 

Вычисляется значение АС, по формуле: 


@+с-—1 
ас,=|| У ас, (@+е-1|+(@Ста +аСшах 2 |/2, (3) 
а 
где АСти и ЯАСтах — минимальный и максимальный элементы вектора 4С, 


соответственно. 
На основе значения С, при просмотре элементов вектора А4С с конца 


выбирается такой индекс „ах, для которого выполняется условие: 

ас; „ 24Са. (4) 
Элементы вектора 4С с индексами от 1 до йлах +1 нормируются так, чтобы 
их сумма была равна 1. 


Формируется изображение / ? по формуле: 


Трах +1 
Е >. ГС ру хе]. (5) 
ЛЕ 
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При формировании цветного результирующего изображения, кроме выполнения 
пунктов а) — к), для каждой цветовой компоненты исходного снимка (т.е. в 


описанных выше формулах 4=1) изображения / ? выполняется следующее 


преобразование: 
Ок Р+1 г (-Р)уЕЬ хехе АЙ, (6) 


причем коэффициент Р для каждой цветовой компоненты вычисляется 
следующим образом: 


пах +1 
Вы вы. За > аб: (7) 
пах + 1=2 


3. Каждая цветовая компонента изображения / . масштабируются на отрезок [0,1], 
после чего применяется метод эквализации гистограммы. Таким образом и 


1 
формируется выходное изображение /°"”', которое может быть цветным или 
полутоновым. 


Экспериментальные результаты 


Представленный алгоритм визуализации был применен при сегментации 
различных полутоновых и цветных медицинских изображений, примерами которых 
являются снимки, представленные на рис. 1 — 3. 

При проведении экспериментов был использован алгоритм нечеткой 
кластеризации $ЕСМ [6, 7]. При его применении динамическое сжатие и обработка 
нечеткой функции принадлежности не использовались. Перед применением карты 
Кохонена число нечетких кластеров увеличивалось в 4 раза, а также осуществлялось 
упорядочение центроидов на основе матрицы Евклидовых расстояний. Выбор 
исходного числа нейронов после обучения карты Кохонена осуществлялся на основе 
максимального соответствия. При кластеризации изображений на рис. 1 аи рис. 2 а 
с=8, а для снимка на рис. За-— с=6. 

На рис. 1 а приведено полутоновое медицинское изображение — томограмма 
головного мозга, выполненная с целью диагностики наличия гематомы, а также 
определения области ее влияния в случае выявления. Однако диагностирование 
гематомы по исходному снимку затруднено из-за ее расположения на 
низкоконтрастном участке изображения (фрагмент, обведенный прямоугольником), 
а область ее влияния вообще визуально  неразличима. Введение 
рентгеноконтрастного вещества (рис. 1 6) хоть и позволяет более четко выделить 
гематому, однако не выявляет области ее влияния. Визуализация результатов 
нечеткой кластеризации с применением метода на основе сингулярного разложения 
(рис. 1 в) не позволяет четко выделить область влияния гематомы. Применение 
предложенного адаптивного алгоритма (рис. | г) позволяет выделить как саму 
гематому, так и область ее распространения при снижении влияния 
неинформативных составляющих (в частности, дефектов пленки). 
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Рис.1. Рентгеновская томограмма головного мозга: а — исходный полутоновой 
снимок (204х201); 6 — результаты введения рентгеноконтрастного вещества; 
визуализация результатов сегментации: в — методом на основе сингулярного 

разложения; г — предложенным алгоритмом 


На рис. 2 а представлена спин-решетчатая Т] релаксация ЯМР (ядерного 
магнитного резонанса) участка головного мозга с наличием злокачественного 
образования. На снимке присутствует низкоконтрастная область (нижняя левая 
часть). Применение предложенного адаптивного алгоритма (рис. 2 в) позволяет 
получить более высокую детализацию для всего изображения в целом по сравнению 
с использованием метода на основе сингулярного разложения (рис. 2 6). 


а б В 
Рис. 2. Спин-решетчатая Т1 релаксация ЯМР участка головного мозга: а — исходный 
снимок (256х256); визуализация результатов сегментации: б — методом на основе 


сингулярного разложения; в — предложенным алгоритмом 


На рис. 3 а, представлены результаты дерматоскопии, проводимой с целью 
визуального диагностирования меланомы. Одним из основных признаков при 
визуальной диагностике, определяющим наличие меланомы, является наличие бело- 
синих структур. В области интереса (обведена прямоугольником) этого изображения 
присутствует низкоконтрастный фрагмент, что затрудняет однозначную визуальную 
диагностику меланомы и области ее распространения. Применение метода 
сравнения с исходными данными [8] для визуализации результатов нечеткой 
кластеризации (рис. 3 6) хотя и позволяет выделить бело-синие структуры (что, к 
сожалению, не видно из-за полутоновой печати), но из-за излишней детализации 
диагностика области распространения меланомы затруднена. Использование 
предложенного алгоритма адаптивной визуализации (рис. 3 в) позволяет как 
выделить саму меланому, так и область ее распространения при отсутствии 
излишней детализации, затрудняющей визуальную диагностику. При этом следует 
отметить, что применение предложенного алгоритма приводит к некоторому 
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искажению цветов по сравнению с исходным изображением. 


Рис. 3. —- Медицинское КОВ изображение: а — исходный снимок (323х215); 
визуализация результатов сегментации: б — на основе сравнения с исходными 
данными; в — предложенным алгоритмом 


Выводы 


Предложенный в данной работе алгоритм адаптивной визуализации резуль- 
татов сегментации на основе сингулярного разложения позволяет в целом повысить 
чувствительность и достоверность сегментации. При этом предложенный алгоритм 
позволяет для цветных входных снимков получать результирующее цветное 
изображение с достаточно высокой степенью цветового соответствия. К недостаткам 
предложенного алгоритма следует отнести некоторое снижение детализации за счет 
использования исходного снимка в процессе визуализации, что в некоторых случаях 
может усложнить визуальную идентификацию объектов интереса. 
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КЕЗОМЕ 
Г.С. Аспте та, А.А. Уесогоу 
Тье Ме®о4 о{ Ве Адарбуе У15па|7айоп о? Ппасез Зеотет(айоп Кези! 5 Вазед оп 
{пе ЭЗшошаг Уаше$ Эесотро$1 оп 

ТЬ$ агисе 15 деусе4 1ю дезсирноп ап4 ехрегитеш а гезеагсВе$ ог фе а Ш@ез оЁ 
Фе тефо@ оЁ Фе адариуе у15иаПтаНоп ог пптасез зестешайоп гези!$ Базе оп Ше 
зшеаг уаиез десотроз1!юоп. ТВе ргорозе4 а]еогИтт аПо\уз епрапсте ог Ве с1еагпез$ 
ап зепу1уЙу оЁ зеотещаноп о\ше 1ю ашютайоп оЁ фе свотсе оЁ Ше п1о$ з1етисапе 
со[ати$ оР Пе [еЁ зшешаг уесюг$ тайлх ап4 пе теготе у фе пива птазе. 

ТБе ргорозе4 ао шса4ез пех заоез: Ёлтху (пеиго-Ёл7ту) саегие; Фе 
зтещаг уа[лез десотроз оп оР ше таблх Фа сощашз уашез оЁ Ше Био лез$/зресвга! 
сотропепЕ оЁ фе па] птазе ап4 е Ра7ху тетабегзЫр РапсНоп; Ше ащютайс спотсе оРШе 
1105 $1отИЙсапе соитпи$ оЁ фе Пей зшошаг уесюг$ тафлх Базе оп ве соеЁРйс1ет оЁ 1$ 
соитпи$ зеШ-езсирИуепез$ (фе сасшаноп оЁ Шезе соеЁйс1ет и5ез Фе пой зшещаг 
уесюгз таблх); йе Юпише оЁ Фе гези по отаузсае ог со]ог птазе Базе оп фе са|1сёайоп 
оРе \уетомеа зит оЁ\е п1о$ з1етИсапЕ соаии$ о Фе Пей зтешаг уесог$ табпх. 

ТБе ехрегипега гезеагсВе$ оЁ 1е ргорозе а1еогит абШИез оп ехатр/е оЁ Фе 
у15иа17те оЁ Фе зеотещайоп гези$ ог Ше отаузсае ап со]ог теса| ппазез уамои$ 
рнуз1са| пафиге аге зво\уп. 
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Л.Г. Ахметшина, А.А. Егоров 
Метод адаптивной визуализации при сегментации изображений на основе 
сингулярного разложения 

Статья посвящена описанию и  экспериментальному исследованию 
возможностей метода адаптивной визуализации результатов сегментации 
изображений на основе применения сингулярного разложения. За счет 
автоматизации отбора наиболее значимых столбцов матрицы левых сингулярных 
векторов и композиции с исходным снимком, предложенный алгоритм позволяет 
повысить наглядность и чувствительность сегментации. 

Предложенный алгоритм состоит из следующих этапов: выполнение нечеткой 
(нейро-фаззи) кластеризации; сингулярное разложение матрицы, включающей 
значения яркостей цветовых/спектральных компонент исходного изображения и 
нечеткой функции принадлежности; автоматический отбор наиболее значимых 
значений матрицы левых сингулярных векторов на основе коэффициентов 
информативности ее столбцов, которые рассчитываются с использованием матрицы 
правых сингулярных векторов; формирование полутонового или цветного 
выходного изображений на основании вычисления взвешенной суммы наиболее 
значимых столбцов матрицы левых сингулярных векторов. 

Приведены результаты экспериментальных исследований возможностей 
предложенного алгоритма на примере визуализации результатов сегментации 
полутоновых и цветных медицинских изображений различной физической природы. 
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